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DIAGNOSIS  OF  SURFACE  WEATHER 


CONDITIONS  FROM  OBSERVED  AND 


PROGNOSTIC  UPPER-AIR  PARAMETERS 


ABSTRACT 


Objective  techniques  are  being  developed  for  interpreting  grid-point  analyses 
and  prognoses  produced  by  computerized  dynamical  models  in  terms  of  concomitant 
surface-weather  conditions.  This  Technical  Report  describes  the  project  and  work 
accomplished  on  it  since  May  1963. 

Multiple  regression  equations  were  derived  to  express  statistical  relation¬ 
ships  between  surface -weather  variables  and  derived  upper-air  parameters  repre¬ 
senting  pertinent  physical  processes  taking  place  between  the  surface  and  the  500-mb 
level.  These  upper-air  (predictor)  parameters  were  derived  from  observed  height 
and  thickness  values  and  the  climatological  statistics  of  these  values. 

The  work  presently  being  conducted  and  plans  for  future  work  are  discussed. 
Improvement  is  being  sought  by  the  definition  of  better  predictor  parameters  to 
represent  orographic  effects  and  by  the  incorporation  of  moisture  (cloud  amount) 
information  now  available  from  dynamical  models.  The  equations  will  be  tested  on 
real-time  upper-air  prognoses  and  readied  for  use  in  an  operational  test  by  the  Air 
Weather  Service  by  September  1965. 
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SECTION  I 


INTRODUCTION 


1.  Purpose 

This  is  an  interim  Technical  Report  on  Project  3.0  of  the  433L  Meteorological 
Technique  Development  Program  (Cohtract  AF  19(628)-3437).  Objective  techniques 
arc  being  developed  to  interpret  grid-mesh  analyses  and  prognoses  produced  by 
computerized  dynamical  models  in  terms  of  concomitant  surface-weather  conditions 
such  as  ceiling,  visibility,  etc.  The  techniques  can  be  applied  to  analyses  to  diag¬ 
nose  current  terminal  weather  conditions  at  locations  for  which  no  surface  observa¬ 
tions  are  available,  and  can  be  applied  to  prognoses  to  obtain  predictions  of  surface- 
weather  conditions  at  selected  terminals.  Previous  work  on  this  project  (reported 
through  May  1963  [2])  indicated  the  feasibility  of  deriving  generalized  statistical 
operators  for  diagnosing  surface-weather  conditions  for  concomitant  upper-air 
parameters.  The  present  report  covers  the  period  from  October  1963  through 
December  1964. 

2.  Objectives 

The  earlier  report  [2]  described  results  obtained  with  statistical  regression- 
equation  models  whose  predictor  terms  were  restricted  to  parameters  derivable 
from  500-,  700-,  and  1000-mb  height  values,  the  long  term  climatology  of  these 
fields  (i.e.,  mean  values  and  standard  deviations),  and  the  three  thickness  layers 
defined  by  them.  The  equations,  derived  for  the  period  October  through  April, 
employ  predictors  of  varying  complexity,  intended  to  represent  pertinent  physical 
processes  taking  place  between  the  surface  and  the  500-mb  level  (e.g.,  thickness  and 
vorticity  advection).  Predictor  parameters  were  chosen  subjectively,  primarily  on 
the  basis  of  meteorological  knowledge  and  experience,  but  were  constrained  to  forms 
which  could  be  incorporated  readily  into  a  routine  operational  forecast  system. 

The  work  described  in  this  report  has  had  the  following  specific  objectives: 

(a)  to  expand  the  generalized  regression-equation  system  described  in  [2] 
so  that  equations  are  available  for  all  seasons  of  the  year. 
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(b)  to  enhance  the  diagnostic  capability  of  the  equations  by  using  additional 
predictors  which  could  take  into  account  (i)  vertical  motion  due  to  orographic 
effects,  (ii)  “coastal  effects”  resulting  from  permanent  moisture  sources  (such  as 
oceans  and  lakes)  which  affect  surface  conditions  at  specified  grid  points,  and  (iii) 
general  radiation  effects  definable  from  time-of-day  and  latitudinal  considerations. 

3.  Scope 

The  surface- weather  variables  (predictands)  considered  here  are  ceiling, 
visibility,  total  cloud  amount,  precipitation,  and  IOC  (integrated  operating  conditions, 
i.e.,  categories  of  combined  ceiling-visibility  conditions,  such  as  below  1500  feet  and 
3  miles,  etc.).  Regression  equations  are  being  derived  to  express  statistical  relation¬ 
ships  between  the  predictand  variable  and  the  derived  (predictor)  parameters,  under 
the  “perfect  prog”  concept,  i.e.,  the  value  of  a  predictand  variable  at  a  grid  point  is 
inferred  from  concomitant  values  of  predictor  parameters  observed  at  or  computed 
for  that  point  alone.  Further,  the  equations  are  generalized  in  that  a  single  regres¬ 
sion  equation  is  applicable  at  any  point  over  an  extensive  geographical  region. 
Generalization  is  accomplished  in  two  ways:  (i)  by  using  predictor  parameters  based 
on  anomalies  from  local  climatic  normals  and  (ii)  by  suppressing  statistical  relation¬ 
ships  pertinent  to  specific  points  only  or  by  devising  generalized  predictor  parameters 
representative  of  such  effects  at  many  locations,  such  as,  a  “coastal  effect”  term.^ 

The  equations  are  derived  from  data  for  16  widely  scattered  points  over  the  eastern 
and  central  United  States.  The  applicability  of  the  equations  to  other  areas  was 

assessed  by  applying  them  to  selected  series  of  200-grid-point  synoptic  maps 
over  the  entire  United  States  and  adjacent  areas.  In  future  work,  the  equations 
derived  from  U.S.  data  will  be  evaluated  on  maps  of  other  areas  of  the  hemisphere. 

4 .  Data 

The  five-year  (April  1955  through  March  1960)  upper-air  data  sample  con¬ 
sisted  of  twice-daily  hemispheric  manually  gridded  synoptic  analyses  of  500- ,  700- , 
and  1000-mb  height  available  on  magnetic  tape.  The  surface  data  consisted  of 


For  definition  and  derivation  of  the  “coastal  effect”  term,  see  Appendix  II, 
parameter  XCLL. 
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twice-daily  (OOZ  and  12Z)  observations  for  the  period  April  1955  through  December 
1958  from  31  U.S.  stations  east  of  the  Rocky  Mountains.  Each  station  was  selected 
to  be  as  near  as  possible  to  an  upper-air  (standard  JNWP)  grid  point.  Data  from  16 
stations  were  combined  to  derive  the  generalized  equations.  Data  from  the  remain¬ 
ing  15  stations  were  reserved  for  testing  the  generalized  equations.  Figure  1  shows 
the  location  of  the  31  surface  stations. 

Climatological  statistics  required  for  the  computation  of  anomaly-based 
predictors  were  derived  from  the  full  5-year  upper-air  data  samples.  A  discussion 
of  the  procedure  used  was  reported  previously  [2].  These  climatological  data  were 
compiled  on  magnetic  tapes,  by  months. 

The  computations  of  all  derived  predictor  parameters  were  performed  for 
each  in  the  data  sample  before  the  statistical  screening  procedure  was  applied. 

Daily  climatological  statistics  were  obtained  by  interpolation  between  monthly  values. 
Each  case  in  the  data  sample  consists  of  concomitant  observations  of  5  surface 
(predictand)  variables  and  81  predictor  parameters  at  one  specific  station  and  Its 
nearby  grid  point.  All  cases  from  the  16  stations  (see  Figure  1)  were  combined 
into  one  generalized  data  sample  to  derive  the  equations  (one  such  sample  was 
compiled  for  each  season).  Similarly,  a  generalized  data  sample  (for  each  season) 
was  compiled  from  the  remaining  15  stations  as  an  independent  test  sample. 


3 


4 


SECTION  II 


SUMMARY 


5.  Background,  Scope,  Results 

Other  techniques  are  being  developed  for  interpreting  grid-point  analyses 
and  prognoses  produced  by  computerized  dynamical  models  in  terms  of  concomit¬ 
ant  surface-weather  conditions.  This  Technical  Report  describes  the  project  and 
work  accomplished  on  it  since  May  1963. 

Multiple  regression  equations  were  derived  to  express  statistical  relation¬ 
ships  between  surface-weather  variables  (such  as  ceiling  and  visibility)  and  derived 
upper-air  parameters  representing  pertinent  physical  processes  taking  place 
between  the  surface  and  the  500~mb  level  (such  as  vorticity  and  thickness  advection). 
These  upper-air  (predictor)  parameters  were  derived  from  observed  height  and 
thickness  values  and  the  climatological  statistics  of  these  values. 

The  equations  are  generalized  in  that  a  single  equation  is  applicable  at  any 
point  over  an  extensive  geographical  region.  The  derived  predictor  parameters 
must  be  recomputed  for  each  specific  point  at  which  an  equation  is  to  be  solved. 

The  generalized  equations  are  discussed  as  statistical  alternatives  to  dynamical 
models  which  may  be  operationally  infeasible  or  as  yet  undeveloped. 

Solution  of  a  given  equation  for  a  specific  point  provides  an  ordered  numerical 
index  which  represents  the  probability  of  occurrence  of  operationally- significant 
categories  of  the  predictand  (surface-weather)  variable.  Probability  distributions 
for  this  purpose  were  obtained  empirically  by  solving  the  equations  for  all  cases 
of  the  developmental  data  sample.  Solutions  of  the  equation  for  many  grid  points 
over  a  synoptic  map  (using  the  independently  derived  predictor  parameters  for 
each  grid  point)  yield  directly  to  a  diagnosis  of  the  concomitant  surface-weather 
patterns. 

Equations  were  derived  separately  for  4  seasonal  periods  for  each  of  5  pre¬ 
dictand  variables:  ceiling,  visibility,  total  cloud  amount,  IOC  (integrated  operating 
conditions;  i.e.,  combined  categories  of  ceiling  and  visibility),  and  the  occurrence 
of  precipitation.  The  equations  were  derived  on  a  sample  of  data  combined  from 
16  surface  stations  over  the  central  and  eastern  United  States  and  tested  on  an 
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independent  sample  from  15  other  stations  in  the  same  region.  The  results  were 
compared  with  the  results  that  would  be  obtained  by  applying  climatological 
frequencies  (probabilities)  in  all  cases.  Improvement  over  the  climatological 
control  technique  ranged  from  8  to  10  percent  for  ceiling,  7  to  14  percent  for  total 
cloud  amount,  3  to  5  percent  for  visibility,  7  to  9  percent  for  IOC,  and  12  to  16 
percent  for  precipitation. 

Two  3-day  series  of  synoptic  maps,  presented,  show  the  application  of  the 
equations  over  a  200-point  grid  network  over  the  entire  United  States  and  adjacent 
areas.  Subjective  examination  of  these  maps  shows  that  the  equations  have  ability 
to;  (1)  indicate  appropriate  changes  in  the  surface-weather  patterns  corresponding 
with  changes  in  the  upper-air  fields,  (2)  indicate,  where  appropriate,  sharp  gradients 
in  the  surface- weather  index  fields,  and  (3)  locate  the  relative  position  of  poor 
surface  conditions  with  respect  to  surface  pressure  systems.  The  equations  are 
shown  to  have  obvious  weaknesses,  in  particular  with  respect  to  poor  surface 
conditions  resulting  from  orographic  effects  over  mountainous  areas. 

Improvement  is  being  sought  by  the  definition  of  better  predictor  parameters 
to  represent  orographic  effects  and  by  the  incorporation  of  moisture  (cloud  amount) 
information  now  available  from  d3niamical  models.  The  equations  will  be  tested  on 
real-time  upper-air  prognoses  and  readied  for  use  in  an  operational  test  by  the 
Air  Weather  Service  by  September  1965. 
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SECTION  III 


APPROACH 


6  Generalized  Statistical  Operators 

The  concept  utilized  in  this  study  is  to  devise  a  single  statistical  operator 
(e  g  a  regression  equation  or  discriminant  function)  which  can  estimate  a  given 
surface-weather  variable  at  any  point  over  an  extensive  geographical  region. 

Most  applications  of  statistical  methods  to  meteorological  prediction  have  involved 
the  analysis  of  historical  data  from  a  single  geographic  point  or  small  area  to 
obtain  a  prediction  system  applicable  to  that  point  or  area  only.  The  derivation 
of  a  generalized  statistical  operator,  on  the  other  hand,  is  approached  through 
the  analysis  of  a  single  data  sample  comprised  of  observations  from  numerous 
points  over  an  extensive  area  In  the  present  work,  for  example,  generalized 
regression  equations  for  diagnosing  surface-weather  conditions  as  a  function  of  other 
parameters  have  been  derived  from  a  combined  data  sample  of  observations  from 
16  stations  widely  scattered  over  the  central  and  eastern  United  States  (see  Figure  1). 
The  basic  justification  for  assuming  that  generalized  statistical  operators  are 
feasible  lies  in  the  concept  that  whereas  the  numerical  values  used  to  describe  a 
given  meteorological  parameter  may  be  strongly  influenced  by  the  climatic  regime 
and  geographical  location  in  which  it  is  observed,  the  same  basic  physical  laws 
apply  over  all  seasons  and  regions  of  the  earth.  The  underlying  approach  in  general¬ 
ized  statistical  analysis  is  to  examine  historical  data  in  such  a  way  as  to  define 
generalized  relationships  which  reflect  the  operation  of  basic  physical  processes 
common  to  all  regions. 

Dynamical  prediction  models  are  the  foremost  example  of  the  use  of  general¬ 
ized  equations  Dynamical  prediction  equations  are,  in  fact,  direct  representations 
of  the  operation  of  known  physical  laws  and  consist  of  mathematical  expressions 
which  simulate  the  processes  in  the  free  atmosphere.  At  present,  however,  there 
are  no  operationally-useful  dynamical  models  available  that  describe  and  predict 
fields  of  critical  surface- weather  parameters  such  as  ceiling,  visibility,  etc.  The 
generalized  statistical  operator  is  conceived  as  a  possible  interim  solution  to  this 
problem  The  types  of  predictor  parameters  examined  (for  the  statistical  equations) 
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are  similar  to  those  which  might  be  used  in  a  true  dynamical  model;  i.e.,  they 
represent  processes  which  meteorological  experience  and  insight  suggest  are 
pertinent  to  the  problem  and  are  capable  of  estimation  from  routinely-observed 
data.  However,  the  choice  of  individual  parameters  for  any  given  equation,  and 
the  weighting  coefficients  to  be  assigned  to  each,  is  determined  by  a  statistical 
analysis  of  historical  data,  rather  than  by  recourse  to  a  preconceived  physical 
model.  The  basic  statistical  nature  of  such  equations  suggests  that  their  solutions 
be  interpreted  in  probabilistic  terms. 

In  the  present  work  generalized  statistical  operators  are  being  derived  as 
part  of  a  mixed  dynamical-statistical  system.  Reliance  is  placed  on  the  free- 
atmosphere  dynamical  model  as  a  predictive  mechanism.  The  statistical  equations 
are  intended  for  use  primarily  as  a  diagnostic  tool;  i.e.,  they  are  derived  to  repre¬ 
sent  concomitant  surface-upper-air  relationships  which  can  be  applied  to  dynamical 
prognostic  charts  and  analyses. 

The  primary  reasons  for  attempting  generalized  statistical  operators  in  lieu 
of  series  of  “point  location”  operators  are  practical  reasons: 

(a)  In  many  areas  of  the  world,  historical  data  samples  are  unavailable  or 
are  too  small  to  obtain  usable  results  by  statistical  analyses. 

(b)  Where  predictions  or  diagnoses  for  many  locations  are  desired,  the 
cost  in  time  and  resources  may  preclude  the  derivation  of  a  multitude  of  specific- 
location  operators. 

The  derivation  of  generalized  statistical  models  for  describing  surface- 
weather  upper-air  relationships  poses  3  basic  problems; 

(a)  Predictor  parameters  must  be  devised  which  are  operationally-useful 
measures  of  the  pertinent  physical  processes.  The  usefulness  of  a  generalized 
statistical  operator  depends  entirely  on  how  well  the  predictors  represent  such 
processes. 

(b)  The  statistical  models  must  have  a  built-in  adaptability  to  varying  climatic 
regimes,  even  though  the  geographical  and  temporal  extent  of  such  regimes  seldom 
can  be  defined  by  simple  mathematical  functions  of  latitude,  longitude,  and  time  of 
year. 


8 


(c)  The  physical  characteristics  of  the  local  area  surrounding  a  given 
observing  station  or  grid  point  (e.g.,  terrain,  presence  of  coastlines,  smoke  sources, 
etc.)  may  give  rise  to  significant  statistical  relationships  which  are  valid  only  for 
that  location.  In  an  objective  statistical  analysis  such  as  the  so-called  “screening” 
procedures,  the  existence  of  such  singular  relationships  may  preclude  the  selection 
of  more  generally  useful  predictor  parameters. 

In  the  present  work  the  3  problems  have  been  attacked  in  the  following  manner: 

(a)  Generalized  predictor  parameters  were  derived  from  the  original  data 
prior  to  the  statistical  analyses.  Predictors  representing  processes  such  as 
vorticity-  and  thickness-advection  are  in  a  form  suitable  for  computation  with 
routinely  available  data.  Where  more  than  one  term  is  feasible  to  represent  a 
given  process,  all  terms  were  computed  and  the  choice  left  up  to  the  objective 
statistical  screening  process. 

(b)  The  elimination  of  climatic  bias  is  attempted  through  the  use  of  anomalies, 
i.e.,  the  deviations  of  height  and  thickness  values  from  long-term  climatological 
normals.  The  consideration  of  anomaly  fields  of  this  type  as  descriptive  and  pre¬ 
dictive  parameters  is  an  established  and  long-used  procedure,  and  is  described 
extensively  in  professional  literature,  e.g.,  [3,  5].  In  the  present  work  the  use  of 
height-  and  thickness-anomaly  fields  has  proved  especially  useful. 

(c)  The  problem  of  geographic,  or  “local”  effects  has  been  approached  in 
two  ways: 

(i)  Generalized  predictor  parameters  were  derived  to  describe 
permanent  local  physical  characteristics  common  to  numerous  locations. 

These  include  terms  to  represent  local  terrain  and  the  resulting  orographic 
effects,  and  to  represent  the  effects  of  nearby  ocean  and  lake  shorelines. 

(ii)  The  application  of  statistical  screening  procedures  to  generalized 
data  samples  (data  from  many  stations)  serves  to  suppress  the  selection  of 
predictor  parameters  pertinent  to  restricted  locations  only. 
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7.  Statistical  Method  Used 

The  Screening  Multiple  Linear  Regression  Technique  [4]  has  been  the  basic 
statistical  tool  in  this  work.  The  technique  is  applied  to  a  data  sample  consisting 
of  observations  of  a  predictand  variable  and  a  relatively  large  set  of  different 
predictor  parameters.  An  efficient  subset  of  predictors  is  selected  objectively 
by  computer  and  a  multiple  regression  equation  is  derived  to  estimate  the  pre¬ 
dictand  variable  as  a  linear  function  of  the  subset  of  selected  predictors.  The 
selection  of  predictors  is  such  that  each  contributes  significant  information  to  the 
estimate  over  and  above  that  contributed  by  all  of  the  other  predictors  selected. 

The  linearity  of  the  equation  resides  only  in  the  weighting  coefficients  assigned  to 
the  selected  predictors.  The  predictors  themselves,  in  particular  those  derived 
to  represent  physical  processes  such  as  vorticity  advection,  may  be  highly  non¬ 
linear  . 

8.  Predictand  and  Predictor  Parameters 

Because  of  the  discontinuities  inherent  in  observations  of  surface-weather 
elements  (e.g.,  unlimited  ceilings  and  visibilities)  the  predictand  data  were  sub¬ 
jected  to  a  normalizing-transformation  procedure  aimed  at  making  such  variables 
more  amenable  to  statistical  regression  analysis.  For  detailed  discussion  of  this 
procedure  see  [2,  Appendix  A].  Precipitation  data  were  available  only  in  dichotomous 
form,  i.e.,  occurrence  or  nonoccurrence  at  the  time  of  observation.  IOC  (see 
section  4)  was  considered  in  4  operational  categories,  i.e.,  (i)  ceiling  below  1500 
feet  and/or  visibility  less  than  3  miles,  (ii)  ceiling  between  1500  and  5000  feet  and/ 
or  visibility  between  3  and  5  miles,  (iii)  ceilings  5000  feet  or  more  (but  not  imlimited) 
and/or  visibility  between  5  and  7  miles,  and  (iv)  ceilings  unlimited  and  visibility 
7  miles  or  greater. 

Numerous  predictor  parameters  were  derived  to  represent  different  atmos¬ 
pheric  processes.  These  included  representation  of  vorticity  concepts,  thermal 
advection  and  changes,  vertical  motion,  and  stability.  All  of  the  parameters  were 
based  on,  and  derivable  from,  observed  500-,  700-,  and  1000-mb  height  values, 
climatological  mean  values  and  standard  deviations  for  those  three  heights,  and 
the  three  thickness  layers  defined  by  those  heights. 
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Note  that,  thus  far,  no  parameters  based  on  direct  observations  of  the  moisture 
field  or  temperature  patterns  have  been  utilized  (thermal  processes  were  approached 
from  thickness  considerations).  At  the  beginning  of  the  project,  data  on  these  ele¬ 
ments  were  not  available  in  usable  form.  Further,  although  the  techniques  are  being 
developed  on  concomitant  surface  and  upper-air  data,  they  are  intended  primarily 
for  use  with  dynamical  upper-air  prognostic  charts  to  infer  future  surface-weather 
conditions.  At  the  start  of  the  work,  dynamical  prognoses  of  upper-air  moisture 
and  temperature  patterns  were  either  unavailable  or  of  questionable  accuracy  for 
the  purposes  of  objective  statistical  analysis. 

More  recently,  however,  the  development  of  dynamical  moisture-  and  temper¬ 
ature-pattern  prediction  models  has  progressed  so  that  the  inclusion  of  such 
parameters  appears  feasible.  Future  plans  for  the  present  work  are  aimed  at  this 
improvement  (see  SECTION  IV). 

A  major  portion  of  the  work  reported  here  involved  the  definition  and  deriva¬ 
tion  of  additional  predictor  parameters  that  can  be  compared  routinely  from  informa¬ 
tion  that  would  be  available  in  an  operational  system.  These  parameters  include: 

(a)  A  parameter  representing  vertical  motion  caused  by  orographic  effects. 
This  parameter  was  computed  by  an  expression  similar  to  advection  computations 
so  as  to  be  proportional  to  the  strength  of  the  geostrophic  flow  at  1000  mb  across 
contours  of  terrain  elevation. 

(b)  A  parameter  to  reflect  so-called  “coastal  effects”;  i.e.,  the  presence  of 
a  nearby  ocean,  lake,  or  other  permanent  moisture  source.  This  parameter  was 
also  defined  by  an  advection-type  expression  which  was  proportional  to  the  com¬ 
ponent  of  the  geostrophic  flow  at  1000  mb  from  the  direction  of  the  designated 
moisture  source  toward  the  grid  point  for  which  surface-weather  conditions  are  to 
be  estimated  (see  Appendix  II,  XCLL). 

(c)  A  number  of  parameters  to  represent  radiation  or  time-of-day  conditions 
at  each  grid  point.  These  expressions  assume  that  the  local  (sun)  time  would  be 
available  in  the  operating  system  for  all  grid  points,  either  in  tables  or  computed 
from  latitude-longitude  considerations.  The  terms  are  relatively  simple;  they 
include  the  length  of  time  between  observation  time  and  sunrise  or  sunset,  the  dura¬ 
tion  of  daylight  or  nighttime,  etc. 
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In  addition  to  the  data  previously  described,  the  parameters  described  above 
require  this  information  in  an  operational  system: 

(a)  Terrain  elevation  values  (preferably  smoothed)  for  each  grid  point 
where  surface  conditions  are  to  be  estimated. 

(b)  An  array  of  “flag”  integers  to  properly  designate  each  grid  point  as  a 
“coastal”  or  “non-coastal”  point.  The  array  of  such  integers  would  be  permanent, 
at  least  through  a  given  seasonal  period. 

(c)  Appropriate  values  of  latitude  and  longitude  for  use  in  converting  observa¬ 
tion  (Z)  time  into  the  local  (sun)  time  at  each  grid  point.  Latitude  and  longitude 
values  could  be  established  as  permanent  arrays  in  computer  storage  or  computed 
when  needed  as  functions  of  the  grid  point  coordinates. 

The  computational  forms  of  all  predictor  parameters  selected  for  the 
regression  equations  are  shown  in  Appendix  II. 
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SECTION  IV 


THE  EQUATIONS 


9.  Derivation 

2 

A  total  of  55  regression  equations  was  derived  to  estimate  the  following 
surface-weather  elements  (predictands); 

Ceiling 

Visibility 

Total  Cloud  Amount 

Precipitation  (occurrence  or  nonoccurrence) 

IOC  (4  categories;  see  section  8) 

Equations  were  derived  separately  for  each  of  4  seasonal  periods: 

January  through  March 
April  through  June 
July  through  September,  and 
October  through  December 

Three  different  equations  were  derived  for  each  predictand  element  in  each  seasonal 
period.  The  first  equation  (Equation  A)  uses  predictors  selected  from  the  original 
set  of  upper-air  parameters  only  (see  [2,  Appendix  IV]).  For  the  second  equation 
(Equation  B)  the  set  of  possible  predictors  to  be  screened  was  expanded  to  include 
the  orographic-  and  coastal-effect  terms.  For  Equation  C  the  set  was  further 
expanded  to  include  the  radiation,  or  time-of-day  parameters. 

The  objectives  in  providing  a  choice  of  3  equations  for  any  given  application 

were: 

(a)  To  provide  flexibility  with  respect  to  the  computer  facilities  required 
for  operational  implementation.  The  complexity  of  the  computer  programs  for 
applying  the  equations  depends  largely  on  the  types  of  predictor  parameters  used. 
The  application  of  the  type  A  equations  would  be  simplest.  The  type  B  and  type  C 


The  precipitation-index  equations  for  the  July  through  September  season  were 
found  to  be  in  error  and  must  be  rederived. 


13 


equations  would  require  successively  more  involved  computer  programs  and 
capabilities  to  derive  the  terrain,  coastal,  and  time-of-day  predictor  (see  section  8). 

(b)  To  evaluate  the  contribution  of  the  terrain,  coastal  and  time-of-day 
parameters  over  and  above  that  obtained  with  the  original  predictors. 

(c)  To  provide  a  ‘^backup”  for  areas  where  the  type  B  and  type  C  equations 
may  be  found  unsuitable.  For  example,  an  equation  containing  terrain-  and  coastal- 
effect  terms  would  be  inefficient  when  used  over  open-ocean  areas  where  such  terms 
are  not  relevant. 

All  regression  equations  are  listed  in  Appendix  I.  Definition  of  the  symbols 
identifying  the  predictor  parameters  and  the  computational  formulas  for  each  are 
given  in  Appendix  II . 

A  total  of  81  possible  predictors  was  screened  for  the  equations.  Of  these, 

55  were  selected  for  use  in  one  or  more  equations.  The  maximum  number  of  pre¬ 
dictors  selected  for  any  one  equation  was  arbitrarily  limited  to  10.  The  smallest 
number  of  significant  predictors  chosen  for  any  one  equation  was  5.  The  frequencies 
with  which  the  different  types  of  predictor  parameters  were  selected  for  the  type  C 
equations  are  given  in  Table  I. 

Certain  features  of  the  manner  of  predictor  selection  are  worthy  of  note: 

(a)  Most  equations  contain  predictors  representing  distinctly  different 
atmospheric  processes.  Further,  in  each  equation,  each  such  process  is  usually 
defined  by  only  one  (at  most,  two)  predictor  parameters.  For  example,  in  most 
type  A  equations  for  ceiling  the  predictors  represent  the  following  parameters: 

vorticity  field 

vorticity  advection 

the  thickness  field 

thermal  (thickness)  advection 

the  height  change  field 

stability  (spring  and  summer  only) 

(b)  The  predictors  representing  orographic  and  coastal  effects  added  significant 
information  to  the  estimation  of  all  predictands  in  one  or  more  seasons. 


14 


TABLE  I 

FREQUENCY  OF  SELECTION  OF  DIFFERENT  TYPES  OF 
PREDICTOR  PARAMETERS  IN  THE  TYPE  C^  EQUATIONS 


Predict  and  Variable 

Type  of 
predietor^ 

Ceiling 

Visibility 

Total 

cloud 

amount 

IOC 

(cig-vsb) 

Precip- 

No 

precip 

Total 

Coastal -effect 

terms 

6 

8 

4 

6 

4 

28 

Time-of-day.  or 
radiation  term 

4 

9 

G 

6 

0 

25 

Space -mean  vorticity 
fields 

G 

3 

4 

4 

5 

22 

Advection  of  thickness 
anomalies 

4 

4 

4 

4 

3 

19 

Advection  of  height 
anomalies 

4 

4 

4 

4 

2 

18 

Terrain-effect 

l^^rm 

4 

2 

4 

4 

3 

17 

Height  change  fields 

3 

4 

4 

4 

1 

16 

Stability  parameters 

4 

0 

2 

3 

3 

10 

Thickness  anomalies 

1 

2 

2 

1 

1 

7 

Height  anomalies 

1 

1 

1 

3 

1 

7 

Height  fields 

1 

1 

1 

1 

3 

7 

Advection  of  thick¬ 
ness  change 

0 

0 

2 

0 

0 

2 

Advection  of  height 
change 

0 

0 

0 

0 

1 

1 

Thickness  field 

0 

0 

1 

0 

0 

1 

Three  t>pe  C  equations  were  available  for  precipitation,  four  for  all  others: 
see  footnote  2. 

^A  list  of  the  individual  predietor  terms  is  given  in  Appendix  IL 
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(c)  The  time-of-day  parameters  added  still  more  significant  information 
to  the  estimation  of  all  predictands  except  precipitation,  and  were  selected  most 
frequently  in  the  equations  for  estimating  visibility. 

(d)  Parameters  representing  vorticity  advection  were  computed  in  two 
ways.  In  the  first  method,  the  relative  vorticity  field  was  defined  by  the  derivation 
of  observed  heights  from  space  mean  height  values.  In  the  second  method,  the 
vorticity  field  was  represented  by  the  height  anomaly  field  (see  [5,  p.  21]).  The 
first  term  is  a  measure  of  the  flow  across  space  mean  contours;  the  second  term 
measures  the  flow  across  climatic  normal  contours.  Mthout  exception,  when  a 
vorticity  advection  predictor  was  selected  for  an  equation  the  term  based  on  the 
height  anomaly  field  was  chosen,  suggesting  that  the  latter  form  may  be  statistically 
more  stable  and  operationally  more  useful.  This  result  supports  similar  findings 
described  in  earlier  work  by  Martin  [5,  p.  54], 

10.  Interpretation 

A  given  equation  solution  is  an  estimate  of  the  parameter  by  which  observed 
categories  of  the  predictand  variable  were  represented  in  the  generalized  develop¬ 
mental  data  sample  (e.g.,  normalizing  transformation  value,  zero-one  parameter, 
etc.;  see  section  8).  This  estimate  (Y)  is  a  numerical  index  which  can  be  inter¬ 
preted  in  terms  of  the  probability  of  occurrence  of  any  given  category  of  the 
predictand,  at  the  point  for  which  the  equation  was  solved.  A  given  value  of  Y 
consequently  provides  an  estimation  of  the  probability  distribution  over  all  pre¬ 
dictand  categories.  A  set  of  probability  distributions  to  relate  different  values  of 
Y  with  a  given  set  of  predictand  categories  can  be  derived  empirically. 

The  index  (Y  values)  can  be  analyzed  directly  on  a  synoptic  map  (see  section 
11)  as  a  subjective  aid  in  defining  present  and/or  future  areas  of  relatively  “good” 
or  “poor”  surface  conditions  over  the  region  to  which  the  generalized  equation  is 
applied.  The  set  of  empirical  probability  distributions  provides  an  objective  means 
for  interpreting  the  numerical  indices  (equation  solutions)  in  terms  of  observable 
values  (categories)  of  the  predictand  variable. 

For  the  present  equations  (derived  on  U.S.  data),  probability  distributions  over 
selected  sets  of  operationally-significant  categories  of  the  predictand  variables 
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were  obtained  empirically  in  the  following  manner.  The  equation  for  a  given 
predictand  was  solved  for  each  case  of  the  developmental  data  sample.  The 
individual  Y  values  were  arranged  in  ascending  order  and  grouped  into  a  manage¬ 
able  number  of  ordered  ranges  of  Y  value.  The  predictand  category  observed  for 
each  individual  Y  value  was  determined  and  a  contingency  table  prepared  to  show 
the  number  of  times  each  predictand  category  was  observed  within  each  Y  range. 
The  relative  frequency  of  a  predictand  category  within  a  given  Y  range  represents 
the  empirical  probability  of  occurrence  of  that  category,  given  that  the  solution  of 
the  equation  falls  within  that  Y  range.  The  spectrum  of  relative  frequencies  for 
that  Y  range  consequently  represents  an  empirical  probability  distribution  over  all 
predictand  categories. 

Contingency  tables  of  this  tjqje  were  compiled  for  each  predictand  variable 
and  season.  For  the  purposes  of  this  report,  tables  were  derived  for  the  type  B 

5 

equations  (see  section  9).  The  stability  of  the  empirical  probability  distributions 
derived  from  the  developmental  data  was  assessed  by  applying  the  equations  to  the 
reserved  sample  consisting  of  data  for  the  15  stations  not  included  in  the  develop¬ 
mental  sample.  Contingency  tables  were  again  compiled  as  described  above. 

A  Brier-Alien  P  score  [1]  was  computed  using  (i)  the  empirical  probability 
distributions  of  each  table  and  (ii)  the  overall  climatological  frequencies  of  each 
table.  In  general,  this  score  (on  the  contingency  tables  compiled  from  the  inde¬ 
pendent  data)  was  better  (lower)  during  the  spring  and  summer  seasons,  but  showed 
the  greatest  improvement  over  climatology  during  the  fall  and  winter  seasons.  The 
improvement  over  climatology  for  the  5  predictand  variables  was  as  follows; 

Ceiling  :  from  7%  (Summer)  to  10%  (Winter) 

Visibility  :  from  3%  (Winter)  to  5%  (Spring) 

Total  cloud  amount  ;  from  7%  (Summer)  to  14%  (Fall) 
Precipitation  ;  from  12%  (Spring)  to  16%  (Fall) 

IOC  :  from  4%  (Summer)  to  9%  (Winter) 

5 

At  the  time  of  this  report,  computer  programs  were  not  yet  available  to  apply 
to  type  C  equations  (with  time-of-day  radiation  parameters)  on  a  grid-point  net¬ 
work  over  wide  areas.  The  tjqje  B  equations  could  be  so  applied,  and  were  applied 
on  a  series  of  synoptic  maps  (see  section  10). 


17 


It  must  be  remembered  that  the  scores  listed  above  represent  results 
obtainable  with  analyses  or  “perfect  prognoses”  only.  The  results  to  be  expected 
when  real-time  prognostic  charts  are  used  have  not  yet  been  determined  (see 
SECTION  V). 

All  of  the  contingency  tables  are  presented  in  Appendix  in.  The  ranges  of 
Y  value  listed  in  the  tables  correspond  to  those  for  which  contours  were  drawn 
on  the  synoptic-map  series  discussed  in  section  11. 

The  equations  derived  from  U.S,  data  are  being  applied  and  their  solutions 
(indices)  analyzed  over  synoptic  maps  of  different  regions  of  the  Northern  Hemi¬ 
sphere,  such  as  the  entire  United  States  and  adjacent  areas  and  areas  of  Europe 
and  Western  Asia  (see  section  11),  It  appears  that  maps  of  this  type  could  provide 
useful  tools  objectively  analyzing  the  large-scale,  synoptic-map  characteristics 
of  surface-weather  patterns  associated  with  the  upper-air  circulation  analyses 
and  prognoses  produced  by  computerized  dynamical  models. 

The  empirical  probability  distributions  derived  from  the  16-station  general¬ 
ized  sample  appear  to  provide  useful  interpretations  of  index  analyses  in  terms 
of  observable  predictand  values  (such  as  ceiling  heights)  throughout  the  areas 
encompassed  by  the  31  stations  to  which  they  have  been  applied,  viz.,  the  central 
and  eastern  United  States.  The  statistical  stability  of  the  equations  can  be 
assessed  by  a  comparison  of  the  relative  frequency  distributions  and  Brier-Alien 
P  scores  obtained  on  the  contingency  tables  for  the  independent  15-station  sample 
with  those  computed  on  the  16-station  developmental  sample  (see  Appendix  HI). 

In  previous  work  [2],  generalized  equations  derived  from  a  combined  data  sample 
from  10  stations  were  applied  to  5  other  stations:  these  results  were  then  compared 
with  results  obtained  by  “single  station”  equations  derived  specifically  for  each  of 
the  5  stations  from  its  own  data  only.  The  comparison  showed  no  consistent  sig¬ 
nificant  difference  in  the  results. 

The  empirical  probability  distributions  derived  from  the  31-station  U.S.  data 
samples  might  also  be  useful  as  “first  estimates”  for  other  regions,  particularly 
regions  for  which  historical  data  are  sparse  or  unavailable.  Nevertheless,  better 
results  could  undoubtedly  be  obtained  from  additional  relative-frequency  contingency 
tables  compiled  for  regions  or  specific  locations  of  particular  operational  importance. 
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The  form  of  such  tables  (i.e.,  the  number  and  limits  of  predictand  categories) 
could  be  chosen  so  as  to  satisfy  best  the  operational  requirements  of  the  particular 
region  or  station.  It  is  suggested  that  a  valuable  objective  procedure  might  be 
devised  for  forecasting  surface  conditions  at  specific  important  terminals  by  com¬ 
piling  a  scatter -diagram  to  relate  computed  index  values  from  the  generalized 
equations  to  concomitant  values  of  the  predictands  at  that  station, 

11.  Application  to  Synoptic  Maps 

The  type  B  equations  for  the  January  through  March  and  July  through 
September  seasons  have  been  applied  to  three  years  of  daily  (12Z)  upper-air 
charts  for  February  (1956—58)  and  August  (1955—57),  Indices  for  the  5  predictand 
variables  were  computed  for  each  day  over  two  200-pomt  grid  networks,  the  first 
over  the  United  States  and  adjacent  areas  (see  Figure  1),  the  other  over  Europe, 
the  Mediterranean,  and  western  Asia,  Two  3-day  series  of  these  maps  over  the 
U,S.  grid  network  (14—16  February  1958  and  12—14  August  1955)  are  presented 
in  this  report.  The  series  were  selected  as  exhibiting  prominent  and  well-defined 
surface-weather  patterns  which  would  provide  a  means  for  a  subjective  evaluation 
of  the  ability  of  the  equations  to  indicate  such  patterns.  The  February  series  shows 
the  progression  of  a  well-defined  low-pressure  center  from  the  Texas  Panhandle 
across  the  Gulf  Coast  States  and  up  the  Atlantic  coast  to  Long  Island.  A  high- 
pressure  cell  intrudes  southward  into  the  Plains  States,  while  a  typical  occluded 
frontal  system  moves  in  from  the  Pacific  Ocean  across  the  Northwestern  States. 

In  the  August  series,  an  elongated  high-pressure  cell  over  the  central  and  north 
central  United  States  moves  rapidly  eastward  as  a  hurricane,  located  over  Cape 
Hatteras  on  August  12,  moves  inland  to  become  an  extratropical  low  center  over 
Lake  Huron  on  August  14. 

The  maps  for  the  February  series  are  presented  in  Appendix  IV;  those  for 
the  August  series,  in  Appendix  V.  The  following  maps  are  shown  for  each  day  of 
each  series; 

(a)  Maps  of  the  500-mb  height  and  surface-pressure  pattern  analyzed  for 
values  at  grid  points  only.  These  two  fields,  in  conjunction  with  the  700-mb  height 
field  and  height-  and  thickness-climatological  data  (not  shown)  represent  the  basic 
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data  from  which  derived  predictor  parameters  were  computed  to  solve  the 
regression  equations. 

(b)  Separate  maps  for  each  of  the  5  predictand  indices,  analyzed  for  the 

6 

values  computed  at  grid  points. 

(c)  The  observed  surface- synoptic  chart  for  that  day  reproduced  from 
the  Historical  Daily  Series  [6]. 

Detailed  discussion  of  each  map  is  impractical  within  the  scope  of  this  report. 
Ceiling  and  visibility  data  are  not  reported  directly  on  the  observed  surface 
chart,  although  indications  of  value  can  be  inferred  from  the  coded  cloud  heights 
and  present  weather  symbols  shown.  Total  cloud  amount  must  be  inferred  from 
the  shaded  station  circles.  Nevertheless,  the  maps  do  permit  a  subjective  evalua¬ 
tion  of  the  information  furnished  by  the  index  fields.  The  contours  drawn  for  each 
index  field  can  be  interpreted  through  the  appropriate  probability  distributions 
obtained  from  the  contingency  tables  given  in  Appendix  III. 

The  encouraging  features  described  in  the  previous  work  [2]  are  again 
evident  on  the  present  maps: 

(a)  The  equations  show  ability  to  relate  changes  in  the  upper-air  pattern 
to  corresponding  changes  in  the  indicated  surface-weather  patterns.  This  can 
be  observed  in  the  February  series,  for  example,  in  the  movement  of  the  low- 
ceiling  indications  to  correspond  with  the  movement  of  the  500-mb  trough  and 
low-pressure  center  from  Texas  east  and  northeastward. 

(b)  There  is  apparent  skill  in  locating  the  position  of  poor  surface  conditions 
relative  to  surface-pressure  systems.  On  the  maps  for  14  February,  for  example, 
the  indications  for  poorer  surface-weather  conditions,  i.e.,  low  ceilings  and 
visibilities,  (higher  precipitation  probabilities)  are  found  on  the  eastern  side  of 
the  low  surface-pressure  center  over  the  Texas  Panhandle. 

(c)  The  equations  have  an  apparent  ability  to  define,  within  relatively  short 
geographical  distances,  sharp  distinctions  between  areas  of  poor  and  good  surface 

0 

Maps  of  the  precipitation  index  for  the  August  series  are  omitted.  See  foot¬ 
note  2. 
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conditions.  This  is  particularly  evident  on  the  maps  for  August  12  in  the  sharp 
gradient  of  the  ceiling  index  in  the  direction  from  which  the  hurricane  is  moving. 
This  gradient  corresponds  well  with  the  conditions  shown  on  the  observed  surface 
map,  and  contrasts  sharply  with  the  relatively  flat  ceiling-index  gradient  extending 
to  the  northeast,  or  leading  direction  of  the  storm. 

The  incorporation  into  the  equations  of  orographic-  and  coastal-effect 
terms  has  improved  the  results,  in  particular  along  the  west  coast  and  eastern 
slopes  of  the  Rocky  Mountains.  The  type  A  equations  applied  in  the  previous 
report  [2]  did  not  provide  indications  of  low  ceilings,  for  example,  along  the  west 
coast  to  correspond  with  conditions  such  as  were  observed  there  on  14  February. 
Nor  could  they  define,  as  does  the  ceiling-index  map  for  15  February  (derived 
with  the  type  B  equation),  indications  of  orographic  low-cloud  conditions  over 
western  Kansas  and  Colorado. 

The  maps  show  the  equations  to  have  obvious  weaknesses,  in  particular 
when  they  are  applied  over  mountainous  regions.  On  the  maps  for  14  February, 
for  example,  a  low  precipitation  index  (low  probability  of  occurrence)  over  western 
Montana  coincides  with  an  extensive  area  of  continuous  snowfall.  Examination  of 
the  individual  predictor  parameters  in  this  case  showed  that  the  estimation  was 
strongly  influenced  by  the  orographic -effect  term.  In  this  instance,  this  param¬ 
eter  indicated  a  spurious  downslope  flow  resulting  from  unrepresentative  elevation 
values  assigned  to  the  nearby  grid  points.  The  plans  for  future  work  to  improve 
this  and  other  observed  weaknesses  are  discussed  in  the  following  section. 
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SECTION  V 


FUTURE  WORK 


The  work  in  this  project  presently  being  conducted  or  planned  for  the 
immediate  future  is  aimed  at  preparing  a  regression-equation  system  which  can 
be  applied  on  an  operational-test  basis  at  the  Global  Weather  Central  (GWC). 

The  plans  provide  for  the  system  to  be  complete  and  ready  for  testing  by  Septem¬ 
ber  1965.  The  specific  objectives  of  the  work  to  be  accomplished  by  that  date  are: 

(a)  To  improve  the  present  equations  (or  derive  new  equations)  to  utilize 
information  available  from  analyses  and  prognoses  of  upper-air  moisture  patterns 
provided  by  the  dynamical  (backward-trajectory)  cloud  model  now  in  operation  at 
the  Global  Weather  Central.  A  simple  approach,  again  using  the  “perfect  prog¬ 
noses”  concept,  is  being  tried  first.  New  equations  for  ceiling,  visibility,  precipita¬ 
tion  and  IOC  are  being  derived  under  the  assumption  that  the  total  cloud  amount 
observed  at  each  grid  point  is  known  or  has  been  accurately  forecast.  In  application, 
this  predictor  would  be  obtained  from  cloud  analyses  (in  some  regions,  perhaps, 
based  on  satellite  observations)  or  from  the  dynamical  cloud  prognoses. 

(b)  To  refine  the  terrain-elevation  parameters  or  to  re-derive  the  oro¬ 
graphic-effect  parameters  to  eliminate  the  inconsistencies  observed  when  the 
present  equations  were  applied  over  mountainous  areas  (see  section  11).  A 
smoothed  terrain-elevation  field  used  for  the  GWC  dynamical  models  is  being 
examined  as  a  possible  improvement.  In  addition,  the  terrain-effect  parameter  is 
being  re-defined  to  utilize,  for  the  advection-type  computations,  a  measure  of  the 
geostrophic  flow  at  the  elevation  of  the  central  grid  point,  rather  than  a  measure 
based  on  1000-mb  height  values  exclusively.  These  changes  in  this  parameter  can 
be  expected  to  require  re-derivation  of  all  equations  in  which  it  is  used. 

(c)  To  transform  the  present  equations,  or  to  derive  new  equations  which 
can  make  use  of  surface  predictand  (or  predictor)  data  available  for  times  pre¬ 
ceding  the  valid  time  of  the  dynamical  prognoses  being  used  to  solve  the  equations. 
One  way  in  which  the  present  equations  can  be  transformed  for  this  purpose  is  as 
follows: 
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Consider  any  one  of  the  present  equations  (assume,  for  simplicity,  that  it 
uses  only  1  predictor).  It  was  derived  in  the  form 


=  a  +  bX„ 
0  0 


(IV-1) 


where  is  the  predictor  value,  the  index  estimate  for  the  same  time  (time  0) 
and  a  and  b  are  weighting  coefficients.  The  equations  are  applied  to  dynamical 
prognostic  charts  in  the  following  sense: 


Y,j,  =  a  +  bX^ 


(IV-2) 


where  X^  is  the  prognostic  value  of  the  predictor  at  some  time  T  (e.g.,  24  hours) 
in  the  future  which  represents  the  valid  time  of  the  prognostic  chart. 

Assume  that  in  a  given  application  we  have  available  observed  or  analyzed 
values  of  Y^  and  X^,  as  well  as  X.^,.  From  Equation  (IV-1)  we  also  have  Y^. 

Then  Y^  -  Y^  represents  the  diagnostic  error  incurred  by  Equation  (IV-1)  at  the 
given  grid  point  at  initial  time  0.  If  we  assume  that  this  diagnostic  error  persists 
in  its  entirety  throughout  the  forecast  period  from  0  to  T,  we  have  a  revised 
estimate  of  Y  for  time  T,  vis.  Y^,  such  that 


Y  =  Y  +  (Y  -  Y  ) 

T  T  ^0  0^ 

Substituting  Equations  (IV-1)  and  (IV-2)  we  have 

a  -  bX, 


(IV-3) 


Y  =  a  +  bX^  Y^ 


0 


=  Y^ 
T  0 


b(X^  -  X^) 

V  T  o' 


(IV-4) 


or 


Y_  -  Y^  =  b  (X^  -  X^) 
T  0  ^  T  o' 


Equation  (IV-4)  represents  a  forecast  of  the  change  in  the  index  from  its  observed 
value  at  time  0,  based  upon  the  change  forecasted  for  the  predictor  X  for  that 
period  by  the  dynamical  model.  The  treatment  of  additional  predictors  would  be 
the  same.  Under  the  assumption  that  the  initial  diagnostic  error  persists,  no 
re-derivation  of  equations  would  be  required,  i.e.,  the  value  of  b  is  the  same  in 
Eqs.  (IV-1)  and  (IV-4).  The  validity  of  the  basic  assumption  is  being  assessed. 
The  present  equations  are  applied  in  the  revised  form  to  the  data  samples 
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previously  used,  using  values  of  and  observed  12,  24,  and  36  hours  previous 
to  the  time  of  the  current  data.  The  basic  assumption  will  obviously  be  question¬ 
able  for  that  portion  of  the  diagnostic  error  (Y^  -  Y^)  due  to  dynamical  effects 
such  as  poor  predictor  definition  or  the  incorrect  movement  of  synoptic  features 
(e.g.,  troughs  and  ridges)  during  the  forecast  period.  Nevertheless,  it  may  prove 
valid  at  grid  points  where  the  initial  diagnostic  error  is  largely  a  result  of  persistent 
local  effects  (e.g.,  terrain  characteristics)  not  considered  by  the  equations.  It  is 
worth  noting  that  Eq.  (IV-4)  could  be  used  in  a  more  general  sense: 

=  Y^  +  h(X^  -  X^)  (IV-5) 

in  which  the  subscript  t  refers  to  any  given  time  (e.g.,  hour)  during  the  forecast 
period  between  0  and  T.  Equation  (IV-5)  would  provide  a  means  for  utilizing 
asynoptic  observations  of  Y^  and/or  X^  to  adjust  the  original  estimate  Y.^,  to  later- 
observed  data.  Values  of  X^  could  presumably  be  obtained,  if  necessary,  from 
hourly  iterations  of  the  dynamical  prognostic  model. 

(d)  To  evaluate  the  usefulness  of  the  equations  when  applied  to  real-time 
prognostic  data..  In  their  present  form,  the  equations  are  an  optimized  estimate 
of  concomitant  relationships  between  predictand  and  predictor  parameters.  As 
such,  they  are  appropriate  means  for  diagnosing  surface-weather  conditions  over 
areas  for  which  upper-air  analyses  are  at  hand  but  surface  observations  are 
unavailable.  Their  suitability  for  application  to  dynamical  prognostic  charts,  how¬ 
ever,  can  be  assessed  only  on  a  real-time  basis,  i.e.,  by  actually  applying  them  to 
prognostic  charts  and  analyzing  the  results  from  both  objective  and  subjective 
standpoints.  The  complexity  of  the  parameters  utilized  in  the  dynamical  model 
and  the  predictors  used  in  the  regression  equations  would  make  a  purely  analytical 
analysis  of  potential  errors  impractical.  One  obvious  source  of  error  would  be 
the  "flattening”  or  "smoothing”  of  wave  patterns  evidenced  on  dynamical  prognoses 
for  longer  periods,  e.g.,  from  48—96  hours.  Height  anomaly  fields,  for  example, 
could  be  expected  to  be  less  well  defined  on  such  maps.  Separate  equations  could 
be  derived,  of  course,  for  various  forecast  periods  from  real  prognostic  data. 

Such  equations  would,  however,  be  "locked  in”  on  the  particular  dynamical  model 
which  produced  the  prognoses,  and  might  require  re-derivation  to  accommodate 
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changes.  A  possible  simpler  solution  can  be  suggested.  In  lieu  of  separate  equa¬ 
tions  for  different  forecast  periods,  separate  contingency  tables  (i.e.,  probability 
distributions— see  Appendix  UI) could  be  compiled  for  each  pertinent  forecast  period, 
so  that  the  initial  set  of  equations  could  be  applied  and  correctly  interpreted  over 
all  periods.  In  the  event  of  changes  in  or  revisions  to  a  given  dynamical  model, 
the  compilation  of  new  tables  for  longer  forecast  periods  would  be  a  clerical  task 
which  could  commence  immediately  upon  initiation  of  the  changed  model. 

A  test  of  the  present  equations  on  a  sample  of  real-time  prognostic  data 
will  be  conducted  before  completion  of  this  project.  Magnetic  tapes  containing 
dynamical  prognostic  data  for  an  extended  period  have  been  obtained  from  the  • 
Global  Weather  Central  for  the  purpose.  Nevertheless,  it  appears  likely  that  for 
practical  (data-processing)  reasons,  a  comprehensive  test  of  this  nature  must 
await  the  implementation  of  the  equations  on  an  operational-test  basis. 
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APPENDIX  I 


EQUATIONS 


Specification  Equations  -  Season  I  (Jan  -  Mar) 


Predictor  Set  A:  Utilizing  only  predictors  derived  from  1CX)0-rab,  7CX)--mb, 
and  500-mb  Heights  and  Climatology. 


CIG  (EQ  A) 


=  -  0.036662  -  0.00029588 
+  0.00831 52 
-  0.044286 


(XA5)  +  0.00034378 
(DH7)  -  0.014255 
(V5) 


(XATS5) 

(AT75) 


VIS  (EQ  A) 


=  -  .0032182  -  0.00031182  (XA5)  +  0.00033712  (XATS5) 
+  0.0071536  (DH7)  -  0.0047047  U5) 

-  0.025755  (V5) 


TCA  (EQ  A) 


=  -  0.037024  +  0.00035089  (XA5)  -  0.00034799  (XATS5) 
-  0.0071313  (DH7)  +  0.19419  (ST75) 

+  0.028334  (V5) 


PGP  (EQ  A) 


=  1 .8598 


0.00033048  (XV5T) 
0.001737  (H7) 

0.0028765  (DH7) 


+  0.01  8I4 
+  0.0017149 

-  0.0027809 


(V7) 

(A5) 

(XSTS7) 


IOC  (EQ  A) 


=  3.8579  “  0.00038344 
+  0.010211 
-  0.056292 


(XA5)  +  0.00045339  (XATS5) 
(DH7)  -  0.01923  (AT75) 
(V5) 


7 

Notation  in  parentheses  beside  numerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations 


Season  I  (Jan  -  Mar) 


7 


Predictor  Set  B;  Utilizing  predictors  derived  from  1000-mb,  700-mb,  and 

i)00~mb  Heists  and  Climatology  plus  Coastal  and  Terrain 
Terms . 


CIG  (EQ  B) 

=  -  0.042053  -  0.00021316  (XA5)  -  0.000070704  (XIL) 

-  0.027964  (VS)  +  0.00025944  (XATS5) 

+  0.0070219  (DH7)  -  0.013230  (AT75) 

-  0.044080  (V5)  +  0.0053553  (XCLL) 


VIS  (EQ  B) 

=  0.13937  -  0.00025122  (XA5)  +  0.010517  (XCLL) 

+  0.00026933  (XATS5)  -  0.22763  (cm) 

+  0.0046850  (DH7)  -  0.12957  (STS5) 

-  0.023709  (V5) 


TCA  (EQ  B) 


=  -  0.047930  +  0.00029149  .(XA5)  +  0.000063576  (XEL) 

-  0.00035005  (XATS5)  -  0.0064892  (DH7) 

+  0.17991  (ST75)  +  0.034030  (V5) 

+  0.000091544  (XDTS5)  -  0.0043348  (XCLL) 


PCP  (EQ  B) 


=  2.2995  +  0.000026005  (XEL) 

-  0.00026611  (XV5T) 
+  0.010628  (VS) 

+  0.013076  (V5) 


-  0.0021777  (H7) 

-  0.0031202  (XCLL) 
+  0.0035409  (AT75) 


IOC  (EQ  B) 


=  3.8493  -  0.00027794  (XA5)  -  0.000086709  (XEL) 

~  0.036621  (VS)  +  0.00034546  (XATS5) 

+  0.0085045  (DH7)  -  0.017856  (AT75) 

-  0.055900  (V5)  +  0.0071629  (XCLL) 


7 

Notation  in  parentheses  beside  numeric-'- 1  coefficients  design*' tes  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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n 

specification  Equations  -  Season  I  (Jan  -  Mar)  ' 

Predictor  Set  C:  Utilizing  predictors  derived  fron  lOOO-mb.  7C)0-nib.  and 

500-mb  Heights,  Climatology,  Coastal  and  Terrain  Terms 
plus  Radiation  Terms. 

CIG  Ceq  c) 

<=  -  0.27674  -  0.00021495  (XA5)  -  0.000071179  (XEL) 

-  0.027643  (VS)  +  0.00026123  (XATS5) 

+  0.0069794  (DH7)  -  0.013512  (AT75) 

-  0.044695  (V5)  +  0.016548  (sTM) 

40.0052102  OCCLL) 

VIS  (EQ  C) 

=  -  0.13930  -  0.00023494  (XA5)  4-  0.0089670  (XCLL) 

4-  0.00025385  (XATS5)  -  0.23746  (CID) 

+  0.028747  (STM)  4-  0.0046824  (DH7) 

-  0.14583  (STS5)  -0.025074  (V5) 

4-  0.0055241  (THEAT)  -  0.000028983  (XEL) 

TCA  (EQ  C) 

=  0.25013  4-  0.00029122  (XA5  )  4-  0.000063980  (XEL) 

-  0.00035558  (XATS5)  -  O.OO6546I  (DH?) 

4-0.17886  (ST75)  -  0.018045  (STM) 

4-  0.034157  (V5)  4-  0.000094015  (XDTS5) 

-  0.0041988  (XCLL)  -  0.014229  (TAST) 

PCP  (EQ  C) 

No  change  from  EQ  B. 

IOC  (EQ  C) 

=  4.6438  -  0.00027604  (XA5)  -  0.000089356  (XEL) 

-  0.037390  (vs)  4-  O.OOO335O8  (XATS5) 

4-  0.0082448  (DH7)  -  0.01 8325  (AT75) 

4-  0.025942  (STM)  -  0.056605  (V5) 

40.0072725  (XCLL)  -  0.17151  (SR) 

"^Notation  in  parentheses  beside  numerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  II 


(Apr  -  Jun) 


7 


Predictor  Set  A:  Utilizing  only  predictors  derived  from  1CK)0-mb,  700-iiib, 

and  500-mb  Heights  eind  Climatology. 


GIG  (EQ  A) 


=  -  0.671  -  0.0042739  (XS5)  -  0.049184  (V7) 

+  0.00088101  (XATS7)  +  O.468I4  (STVT) 

+  0.0065718  (DH5)  -  0.00064679  (XA7) 


VIS  (EQ  A) 


=  -  0.026439  -  O.OOO3O83O 
+  0.0065696 
-  0.00018153 


(XA7)  +  0.00074589  (XATS7) 
(DH7)  -  0.0011623  (XAS) 
(XA5)  +  0.13233  (STVT) 


TCA  (EQ  A) 


=  0.65645  +  0.00059721  (XA5)  -  0.00070210  (XATS5) 

+  0.029627  (V7)  -  0.43782  (STVT) 

-  0.0069458  (DH5)  +  0.00017726  (XDTS5) 


PCP  (EQ  A) 


=  1.3404  +  0.018426  (V7) 

-  0.00043517  (XV7T) 
+  0.00045185  (XAS) 


+  0.000029095  (XDH5) 

-  0.0012251  (H7) 

-  0.0000861 41  (XSTBA) 


IOC  (EQ  A) 


“  2.7936  -  0.00043406  (XA5)  +  0.77237  (STVT) 

+  0.0010234  (XATS7)  +  0.0082896  (DH5) 

-  0.049855  (VS)  -  0.0017291  (XAS) 


7 

Notj?tion  in  parentheses  beside  numerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  II  (Apr  -  Jun) 


7 


Predictor  Set  B;  Utilizing  predictors  derived  from  1000-mb,  700-mb,  and 

5CX)--mb  Heights  and  Climatology  plus  Coastal  and  Terrain 
Terms. 


CIG  (EQ  B) 


=  0.56055-0.0037313  (XS5)  -0.032143  (V7) 

-  0.000057752  (XEL)  +  O.OOO637IO  (XATS7) 
+  0.44411  (STVT)  -  0.037991  (VS) 

+0.01 681 5  (XCLL)  -  0.017040  (XCLN) 

+  0.0049610  (DH5)  -  0.00050894  (XA7) 


VIS  (EQ  B) 


=  0.24634-0.35623  (cm)  +0.011532  (XCLL) 

-  0.00057291  (XA7)  +  0.00058206  (XATS7) 
+  0.0048199  (DH7)  -  0.031312  (VS) 

+  0.14370  (STVT) 


TCA  (EQ  B) 

=  1.1688  +  0.00054384  (XA5)  -  0.00057325  (XATS5) 

+  0.0072458  (V7)  +  0.000057912  (XEL) 

-  0.72506  (STVT)  -  0.0058727  (DH5) 

+0.037829  (VS)  +0.00013168  (XDTS5) 

-  0.12644  (XCLD)  +  0.10078  (STBS) 


PCP  (EQ  B) 


=  1  .1614 


+  0.015487 

(V7) 

+ 

0.00001 7775 

(XDH5) 

-  0.0057538 

(XCLL) 

+ 

0.0067115 

(XCLN) 

-  0.00010321 

(XV7T) 

- 

0.001 0696 

(H7) 

+  0.000011389 

(XEL) 

+ 

0.000093397 

(XA5) 

-  0.00011635 

(XATS5) 

- 

O.OOO891 26 

(DH5) 

IOC  (EQ  B) 

=  3.0361 


-  0.00015628  (XA5) 

+  O.OOO8424O  (XATS7) 

-  0.000067687  (XEL) 

+  0.020889  (XCLL) 

-  0.026646  (V5) 


+  0.67465  (STVT) 
+  0.0056702  (DH5) 

-  0.056055  (VS) 

-  0.020162  (XCLN) 

-  0.00065771  (XA7) 


7 

Notation  in  parentheses  beside  nvunerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season 


II  (Apr  -  Jun)'^ 


Predictor  Set  C;  Utilizing  predictors  derived  from  lOOO-mb,  700-mb,  and 

500-mb  Heights,  Climatology,  Coastal  and  Terrain  Terms 
plus  Radiation  Terms. 


CIG  (LQ  C) 


=  -  0.33751  -  0.0087223  (XS5)  -  0.040620  (V7) 

-  0.000059424  (XEL)  +  0.000501 61  (XATS7) 

+  0.35985  (STVT)  -  0.043265  (VS) 


+  O.OI9OI6 
-  0.038529 

VIS  (£Q  C) 

.  =  0.20117  -  0.36784 

-  0.00057230 
+  0.0086705 
+  0.048118 


TCA  (EO  C) 

=  4.0806  +  0.00057416 

+  0.030617 
+  0.23121 

-  0.0044024 

-  0.10552 


PCP  (EQ  C) 


No  change  from  EQ  B. 


IOC  (EQ  C) 


«  3.0278  -  0.00037024 
+  0.00083898 

-  0.000068391 
+  0.0060027 

-  0.020476 


(XCLL)  -  0.018077  (XCLN) 
(HTGV) 


(cm)  +0.011172  (XCLL) 
(XA7)  +  0.00060590  (XATS7) 

(UiEAT)  +  0.0053138  (DH7) 

(RADN)  -  0.026791  (VS) 


(XA5)  -  0.00059412  (XATS5) 
(V7)  +  0.000060034  (XEL) 

(TOD)  -  0.0075543  (DH7) 

(TS7)  +  0.00015211  (XDTS5) 
(XCLD) 


(XA5)  +  0.60124  (STVT) 

(XATS7)  +  0.0054088  (DH5) 
(XEL)  -  0.04812/+  (VS) 
(THEAT)  +  0.022435  (XCLL) 
(XCLN)  +  0.085769  (SS) 


7 

Notation  in  parentheses  beside  numerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  III  (Jul  -  Sep) 


7 


Predictor  Set  As  Utilizing  only  predictors  derived  from  lOOO-mb,  700-rab, 

and  500-mb  Heights  and  Climatology. 


CIG  (EQ  A) 


=  -  0.92106  -  0.0079627 
-  0.036441 
+  0.010391 


(XS5)  +  0.67919 

(VS)  +  0.0010976 
(DH5)  +  0.00074017 


(STVT) 

(XV7T) 

(XATS7) 


VIS  (EQ  A) 


=  27.445  -  0.00048849  (XA7)  +  2.9202  (STVT) 
+0.0085821  (DHS)  -  0.031177  (H7) 

-  0.075517  (STBA)  +  0.028150  (A7) 


TCA  (EQ  A) 

=  0.74302  +  0.0097542  (XS5)  -  0.001 5206  ijXVTT) 

-  0.57700  (STVT)  -  0.011089  (DH5) 

-  0.00095481  (XATS7)  +  0.037742  (VS) 


IOC  (EQ  A) 


=  1.7829  -  0.0094278 

-  0.032687 

+  0.011394 

-  0.021232 


(XS5)  +1.4079  (STVT) 
(VS)  +  0.0012922  (XV7T) 
(DH5)  +  0.00092896  (XATS?) 
(STBA) 


7 

Notation  in  parentheses  beside  nvunerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  III 


(Jul  -  Sep) 


7 


Predictor  Set  B;  Utilizing  predictors  derived  froa  1000-mb,  700-mb,  and 

500-inb  Heights  and  Climatology  plus  Coastal  and  Terrain 
Terms . 


CIG  (EO  B) 


=  -  1  .6395  -  0.007U38  (XS3)  +  1  .0836  (STVT) 

-  0.039847  (VS)  -  -  0.000048402  (XEL) 

+  0.0010225  (XV7T)  +  0.0082219  (DH5) 

+  0.00061782  (XATS7)  -  0.017654  (STBA) 

+  0.0075814  (XCLL)  -  0.026531  (V5) 


VIS  (EQ  B) 


=  11  .105  -  0.36329 

+  0.00094149 

-  O0O35615 

-  0.010322 

+  0.079610 


TCA  (EQ  B) 


=  0.53153  +  0.0091105 
-  0.52426 
-  0.0011200 
+  0.046827 
-  0.13952 


IOC  (EQ  B) 


“  1 .9811  -  0.012624 

-  0,045441 

-  0.021578 

+  0.00085730 
+  0.0045204 


(cm)  -  0.00071604  (XA7) 

(XATS7)  +  0.0075152  (DH7) 

(HS)  +  0.028535  (AS) 

(T75)  -  0.0010535  (EL) 

(S5)  -  0.00038399  (XDTS7) 


(XS5)  -  0.0013463  (XV7T) 

(STVT)  -  0.010440  (DH5) 

(XATS7)  +  0.000053764  (XEL) 

(VS)  +  0.22987  (cm) 

(XCLD)  +  0.00035315  (XDTS7) 


(XS5)  +  1 .3584  (STVT) 

(VS)  +  0.010503  (XCLL) 

(STBA)  +0.011597  (DH5) 

(XATS7)  -  0.15947  (cm) 

(XST75)  -  0.000044505  (XEL) 


7 

Notation  in  parentheses  beside  numerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  III  (Jul  -  Sep) 


Predictor  Set  C;  Utilizing  predictors  derived  from  lOOO-mb,  7CX)-Eib,  and 

500-mb  Heights,  Climatology,  Coastal  and  Terrain  Terms 
plus  Radiation  Terns. 

CIG  (EQ  C) 

=  -  1.6242  -  0.0072761  (XS5)  +  1  .1326  (STVT) 

-  0.041836  (VS)  .  -  0.000054845  (XEL) 

+  0.00088825  (XV7T)  +  0.010035  (DH5) 

+  0.00067114  (XATS7)  -  0.015695  (TEST) 

-  0.020689  (STBA)  +  0.0075324  (XCLL) 

VIS  (EQ  C) 


*  1 .8203  -  0.31750 

(cm)  - 

0,051840 

(TEST) 

+  0.10274 

(PJUDN)  - 

0.00052811 

(XA7) 

+  O.OOO65035  (XATS7)  + 

to 

8 

• 

0 

(DH7) 

+  0.0075244 

(XCLL)  - 

0.033964 

(ATS5) 

-  0.030350 

(HS)  + 

0.029484 

(A5) 

TCA  (EQ  C) 


“  1 .0834  +  0.0093423 

(XS5)  -  0.0014892 

(XV7T) 

-  0.49915 

(STVT)  -  0,010506 

(DH5) 

-  0.00089385 

(XATS7)  +  0.19626 

(TCD) 

+  0.000068788 

(XEL)  +  0.046942 

(VS) 

+  0.17601 

(CID)  -  0.14144 

(STM) 

IOC  (EQ  C) 


=  2.3468  -  O.OO98O3I 

(XS5) 

+ 

1 .7786 

(STVT) 

-  0.035453 

(TEST) 

- 

0.041010 

(STBA) 

+  0.013706 

(XCLL) 

0.094354 

(ss) 

+  0.0010338 

(XATS7) 

+ 

0.0099428 

(DH5) 

-  O.OOOO58O5I 

(XEL) 

- 

0.074030 

(DURD) 

7 

Notation  in  parentheses  beside  numerical  coefficients  designates  select¬ 
ed  predictor.  Ih:edictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  IV  (Oct  -  Dec) 


Predictor  Set  A:  Utilizing  only  predictors  derived  from  1000-mb,  700-nib, 

and  500-mb  Heights  and  Climatology. 


GIG  (£w  a) 


=  -  4.0349  -  0.0093836  (XS5)  -  0.055221  (V7) 

+  0.0089440  USTS5)  +  0.0070299  (DH7) 
-  0.0048697  U5)  .  +  0.0022319  (H5) 


VIS  (EQ  A) 


*=  2.5237  -  0.00040858  (XA7)  +  0.00037646  (XATS7) 
+  0.0056249  (DH7)  -  0.0023697  (H7) 

+  0.00033404  (XV5T)  -  0.018470  (V7) 


TCA  (EQ  A) 

=  6.0394  +  O.OO9906O  (XS5)  -  O.OO6223O  (H7) 

-  0.00028338  (XATS5)  -  0.0082682  (DH7) 

+  0.0093479  U7)  +  0.042096  (V7) 

+  0.000056149  (XDH5) 


PCP  (EQ  A) 

=  1.5044  +  0.012940  (V7)  -  0.00013871  (XV5T) 

-  0.0013941  (H7)  +  0.00018445  (XA7) 

-  0.00016112  (XATS7)  -  0.00037877  (XAS) 

-  0.0011644  (DH5)  +  0.0014433  Us) 


IOC  (EQ  A) 


=  -  0.37232  -  0.00034889  (XA5)  -  0.064971  (V7) 

+  0.00038687  (XATS5)  +  0.0067859  (DH5) 
-  0.0068368  U5)  +  0.0043995  (H7) 


7 

Notation  in  parentheses  beside  numerical  coefficients  design/^tes  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  IV  (Oct  -  Dec) 


7 


Predictor  Set  B:  Utilizing  predictors  derived  from  1CXX)-nib,  700-mb,  and 

500-mb  Heights  and  Climatology  plus  Coastal  and  Terrain 
Terms . 


CIG  (EQ  B) 


-  6.3925  -  0.0066930  (XS5) 

-  0.034411  (VS) 


-  0.000077249  (XEL) 


+  0.0065636 


0.0094689  (XCLL)  -  0.0058720 
+  0.0067977  (DH7)  +  0.0067698 
-  0.026928  (V5)  -  0.16095 


(H7) 
(A5) 
(XSTS5) 
(CID) 


VIS  (EQ  B) 


0.37861  -  0.00033570 

+  0.0091678 

-  0.000036143 

-  0.023610 


(XA7)  -  0.37294  (CID) 

(XCLL)  +  0.0059548  (DHS) 
(XEL)  +  0.00032472  (XATS7) 
(VS)  +0.0041583  (DTS7) 


TCA  (EQ  B) 


=  1.0871  +  0.0085472  (XS5)  +  0.000077067  (XEL) 

-  0.0093168  (H7)  +  0.0012239  UTS5) 

-  0.0097330  (XCLL)  -  0.00023202  (XATS5) 

-  0.0082604  (DH7)  +  0.0080568  U7) 

+  0.037005  (VS)  +0.12514  .  (CID) 


PCP  (EQ  B) 


=  2.1247  +  0.014139  (V7)  -  0.0064270  (XCLL) 

-  0.0020252  (H7)  -  0.00012222  (XV5T) 

+  0.000015770  (XEL)  +  0.0053232  (XCLN) 

+  0.0077387  (VS)  +0.00097715  U7) 

+  0.000077037  (XA7)  -  O.OOOO63OI9  (XATS7) 


IOC  (EQ  B) 


=  -  3.8080  -  0.00022968  (XA5)  -  0.000094660 

-  0.041337  (VS)  +  0.0079807 

+  0.011499  (XCLL)  -  0.0077854 

•-  0.28334  (CID)  +  0.0081017 

+  0.00024515  (XATS5)  -  0.034690 


(XEL) 

(H7) 

(A5) 

(DH7) 

(V5) 


n 

'Notation  in  parentheses  beside  numerical  coefficients  design-tes  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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Specification  Equations  -  Season  IV  (Oct  -  Dec)"^ 


Predictor  Set  C;  Utilizing  predictors  derived  from  1000-mb,  700-inb,  and 

500-mb  Heights,  Climatology,  Coastal  and  Terrain  Terms 
plus  Radiation  Terms. 


CIG  (£Q  C) 


=  -  7.1i)32  -  0.0064206  (XS5) 

-  0.038407  (VS) 

+  0.0094862  (XCLL) 
+0.015615  (STM) 

+  0.0063645  (XSTS5) 


-  0.000079566  (XEL) 

+  0.0071142  (H7) 

-  0.0042809  (A5) 

+  0.0072585  (DH7) 

-  0.17696  (CID) 


VIS  (EQ  C) 


.020895  -  0.00029877 

(XA7) 

-  0.34619 

(cm) 

+  0.0092255 

(XCLL) 

+  0.024241 

(STM) 

+  0.0044247 

(DHS) 

-  0.000041775 

(XEL) 

+  0.00022887 

(XATS7) 

+  6.0048805 

(THEAT) 

-  0.026080 

(VS) 

+  0.16526 

(TBSR) 

TCA  (EQ  C) 


=  9.3657  +  0.0083372 

(XS5) 

+ 

0.000080799 

(XEL) 

-  0.0092903 

(H7) 

0.001 6544 

(ATS5) 

-  0.0090886 

(XCLL) 

- 

0.015882 

(STM) 

-  0.00022846 

(XATS5) 

- 

0.0082772 

(DH7) 

+  0.0076757 

(A7) 

+ 

0.035016 

(VS) 

PCP  (EQ  C) 


No  change  fTOm  Equation  B. 


10c  (EQ  C) 


=  -  4.9455  -  0.00021588  (XA5)  -  0.000098113  (XEL) 

-  O.O468O8  (VS)  +  0.026352  (STM) 

+  0.0087646  (H7)  -  0.0058477  (A5) 

+0.011554  (XCLL)  -  0.30932  (CID) 

+  0.0087292  (DH7)  +  0.00022730  (XATS5) 


7 

Notation  in  parentheses  beside  numerical  coefficients  designates  select¬ 
ed  predictor.  Predictor  notation  is  identified  in  Appendix  II. 
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APPENDIX  II 


NOTATIONS  AND  DEFINITIONS  OF  SELECTED  PREDICTORS 


PREDICTOR 

NOTATION 

DEFINITION  OF  PREDICTOR 

STM 

Local  time;  sun  at  zenith  at  1200. 

STM  -  (Time  in  EST)  -  [(Longitude  -  75'’)/l5  ] 

SR 

Time  of  sunrise  in  local  (sun)  as  for  STM. 

(SR  can  be  obtained  from  tables  and/or  formulae 
available  in  nautical  almanacs) 

ST 

Time  of  sunset:  similar  to  SR.  (ST  was  not  selected  as  a  pre¬ 
dictor.  It  is  listed  to  simplify  definitions  to 
follow.) 

TOD 

TOD  =  1,  if  SR  2  STM  s  ST 

TOD  =  0,  if  SR  >  STM  >  ST 

TEST 

Time  until  sunset;  TEST  =-  ST  -  STM,  if  TOD  =  1 

TEST  =0  ,  if  TOD  -  0 

TAST 

Time  since  sunset;  TAST  -  STM  +  24  -  ST,  if  STM  <  SR 

TAST  0  ,  if  TOD  =  1 

TAST  -  STM  -  ST  ,  if  STM  >  ST 

TBSR 

Time  until  sunrise;  TESR  -  SR  -  STM,  if  STM  <  SR 

TESR=  0  ,  if  TOD-  1 

DURD 

Duration  of  daylight;  DURD  --  ST  -  SR 

HTGV 

Heating  parameter  HTGV  =0  .if  TOD  =  0 

HTGV  -  TASR,  if  TOD  =  1  and  STM  ^  12 
HTGV  -  TEST,  if  TOD  -  1  and  STM  >  12 

RADN 

Radiation  parameter;  RADN  --  HTGV,  if  TOD  -  1 

RADN  -CLGV,  if  TOD  -  0 
where  CLGV  -  TAST  if  STM  >  ST 
=  TESR  if  STM  <  SR 

TREAT 

Total  heating  parameter : 

TREAT  -  0  if  TOD  -  0 

TREAT  =  1/2  (TASR)^  if  TOD  -  1  and 

STM  ^  12 

HS 

TREAT  (12  -  SR)^  -  1/2  (ST  -  STM)^ 

if  TOD  =  1  and 
STM  >  12 

Q 

1000-mb  height 

In  dec  a  feet. 
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H7 

O 

700~mb  height 

H5 

SOO-mb  height^ 

DHS 

g 

24-hour  change  in  1000-mb  height 

DH7 

g 

24-hovtr  change  in  700-mb  height 

DHS 

g 

24-hovtr  change  in  SOO-mb  height 

TS7 

g 

1000— 700-mb  thickness 

T75 

g 

700— SOO-mb  thickness 

DTS7 

g 

24-hovtr  change  in  1000— 700-mb  thickness 

AS 

lOOO-mb  height  anomaly®  AS  =  HS  -  HS® 

A7 

700-mb  height  anomaly®  A7  =  H7  -  H7  ® 

AS 

8  — 9 

SOO-mb  height  anomaly  AS  ~  HS  -  HS 

SS 

1000-mb  height  standardized  anomaly  SS  =  AS/c  (HS) 

SS 

SOO-mb  height  standardized  anomaly  SS  =  AS/cr  (HS) 

ATSS 

1000— SOO-mb  thickness  anomaly  ATSS  -  TSS  -  TSS 

=  AS  -  AS 

AT7S 

700— SOO-mb  thickness  anomaly  AT7S  -  T7S  -  T7S 

=  AS  -  A7 

STSS 

1000— SOO-mb  thickness  standardized  anomaly 

STSS  =  ATSS/a  (TSS)^® 

ST7S 

700— SOO-mb  thickness  standardized  anomaly 

ST7S  =  AT7S/cr  (T7S) 

STBA 

Stability  index  derived  from  anomalies 

STBA  ==  ATS7  -  AT7S 

STBS 

Stability  index  derived  from  standardized  anomalies 

10 

STBS  =  STS7  -  ST7S 

8. 

In  deca  feet. 

9.  The  mean  and  standard  deviation  of  a  variable  (at  a  given  time  and  place)  are 
indicated  respectively  by  a  bar  and  the  symbol  a.  e.g.,  HS,  cr  (HS). 

10. 

TSS  represents  thickness  1000  mb  to  SOO  mb 

ATS7  represents  thickness  anomaly  1000  mb  to  700  mb 

ATSS  represents  thickness  anomaly  1000  mb  to  500  mb 

STS7  represents  thickness  standardized  anomaly  1000  mb  to  700  mb. 

40 

STVT 

Stability  index  derived  from  assumed  mean  vertical  temperatures 

STVT  =  0.02866TS7  -  0.03038T75 

VS 

Relative  vorticity  1000  mb^^ 

V7 

11 

Relative  vorticity  700  mb 

V5 

Relative  vorticity  500  mb^^ 

XDH5 

12 

Advection  of  DH5  in  the  500-mb  flow 

XV7T 

12 

Advection  of  V7  in  the  1000— 500-mb  thickness  field 

XV5T 

12 

Advection  of  V5  in  the  1000— 500~mb  thickness  field 

XDTS5 

Advection  of  DTS5  (i.e.,  24-hour  change  in  1000— 500-mb  thick¬ 
ness)  in  the  mean  1000— 500-mb  flow^^ 

XDTS7 

Advection  of  DTS7  (i.e.,  24-hour  change  in  1000— 700-mb  thick- 

12 

ness)  in  the  mean  1000— 700-mb  flow 

XAS 

12 

Advection  of  AS  (i.e.,  1000-mb  height  anomaly)  in  the  1000-mb  flow 

XA7 

12 

Advection  of  A7  (i.e.,  700-mb  height  anomaly)  in  the  700-mb  flow 

XAS 

12 

Advection  of  A5  (i.e.,  500-mb  height  anomaly)  in  the  500-mb  flow 

XS5 

Advection  of  S5  (i.e.,  500-  mb  height  standardized  anomaly)  in  the 

12 

500-mb  flow 

XATS5 

Advection  of  ATS5  (i.e.,  1000— 500-mb  thickness  anomaly)  in  the 

12 

mean  1000— 500-mb  flow 

XATS7 

Advection  of  ATS7  (i.e.,  1000— 700-mb  thickness  anomaly)  in  the 

mean  1000— 700-mb  flow 

11.  The  relative  vorticity  at  a  point  was  computed  by  subtracting  the  space-mean 
value  (obtained  from  the  height  values  at  the  4  surrounding  JNWP  grid  points)  from 
the  height  value  at  the  given  point. 

12.  The  advection  of  a  variable  X  is  computed  from  4  surrounding  points  A,  B, 

C,  D  as  follows  (where  HA  is  the  height  at  point  A) 

Advection  of  X  [  H(A)  -  H(B)]  [X(C)  -  X(D)]  -  [H(C)  -  H(D)] 

[X(A)  -  X(B)] 
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XSTS5 

Advection  of  STS5  (i.e.,  1000— 500-mb  thickness  standardized  anomaly) 

12 

in  the  mean  1000— 500-mb  flow 

XSTS7 

Advection  of  STS7  (i.e.,  1000— 700-mb  thickness  standardized  anomaly) 

12 

in  the  mean  1000— 700-mb  flow 

XST75 

Advection  of  ST75  (i.e.,  700— 500-mb  thickness  standardized  anomaly) 

12 

in  the  mean  700— 500-mb  flow 

XSTBA 

Stability  change  from  advection  of  thickness  anomalies 

XSTBA  =  XATS7  -  XAT75 

XEL 

Index  of  orographically  induced  vertical  motion.  Computed  from 

12 

the  advection  formula  ,  where  X  represents  the  terrain 
heights  at  the  4  surrounding  points. 

CID 

A  zero-one,  or  dummy,  variable  indicating  that  the  grid  point  is 
or  is  not  influenced  by  a  major  moisture  source.  CID  =  1 
if  the  point  is  influenced  by  a  major  moisture  source.  CID  =  0 
if  the  point  is  not  so  influenced. 

XCLL 

Coastal  effect  linear : 

If  CID  =  0,  SCLL  =  0 

If  CID  =  1,  the  4  surrounding  grid  points  are  assigned  values 
of  1  or  0  according  to  whether  there  is  or  is  not  a  moisture 
source  (lake,  ocean,  etc.)  in  the  direction  of  that  grid  point 
which  effects  the  central  grid  point.  XCLL  is  then  computed 
as  the  “Advection”  of  the  assigned  values  by  the  1000-mb 
flow.  If  XCLL  <  0,  the  1000-mb  flow  has  a  component  from 
the  moisture  source  toward  the  central  grid  point.  If  XCLL  >  0, 
the  component  is  away  from  the  central  grid  point. 

E.g., 

0  Or  2 -  coast  line 

•  ^ 

D  / 

0  /  1 

•  oY  • 

k  V  B 

• 

c 

XCLL  =  [(HS^  -  HSg(  (1  -  0)]  -  [(HS^  -  HS^^)  (0  -  1)] 

The  advection  of  a  variable  X  is  computed  from  4  surrounding  points  A,  B, 
C,  D  as  follows  (where  HA  is  the  height  at  point  A) 

Advection  of  X  =  [  H(A)  -  H(B)]  [X(C)  -  X(D)]  -[H(C)  -  H(D)] 

[X(A)  -  X(B)] 
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XCLN 


XCLD 


Coastal  effect,  segmented;  If  XCLL  ^  0,  XCLN  =  XCLL 

If  XCLL  <  0,  XCLN  =  0 

A  zero-one,  or  dummy,  variable  indicating  on-shore  flow  (0)  or 
off-shore  flow  (1). 


43 


TABLE  n 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  CEILING 


APPENDIX  III 


CONTINGENCY  TABLES 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better, 
scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental- sample  table. 
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CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  CEIUNG 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

F  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  VISIBILITY 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better 
F  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  tabic. 
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P  seorc  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

F  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

F  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental- sample  table. 
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< 


c 

"O 

c 

0^ 

a 

o 

TS 

B 


I 

a 


m 

0) 

cd  > 
o  o 

3 

o  ^ 


H  ^ 
d)  U 
>  ^ 
.1^  01 

1 1 


O  d, 

•  w 
o  cd 
Z  o 


"s:  a 

^  I 


o 

Z 


P 

|l 

(2 


o 

z 


o  S 
d  S 

Z  W 


o  B 

r. 


s 


CO 

00 


o 


<£> 


■§ 

O 

U 

a 


o 

ed 

i 


h£> 

B 


I  a 


CO 

m 


i  .s 


Q 

> 

I 


I  a 


s  I 

p 

03  > 

03  ed 
®  £ 

O 

03 

^  .ti 

o  ::: 

<  I 

.  o 
•  a« 

"o  .s  . 

-  w  . 


3 


03 


3 

a  a 

^  I 

2  ^ 

^  a  _ 

— ^  ’ts  d  B  ° 

03  O  = 

03  d3 

^  I 

o  3 

o  X 

i  I 


eif 

cd  o 

a  2  -2 

r  & 

o  ™  4) 

O  ii  E 

a  g. 


03 
> 

iS 

<d  h£) 

TJ  B 


< 

O 

H 


|i 

.M  dJ 

4,-§ 

b£ 

I  ^ 

M 

§  a 

"  n 

03  03 

u  u 
o  o 


3 

•S 

_  .  a 

8  ii  -g 
10,  |0< 


50 


TABLE  Vni 

CONTINGENCY  TABLE  OF  INDEX  VALUES  vs.  OBSERVED  VALUES  OF  IOC 
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TABLE  IX 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  IOC 
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Q 


I 
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Total  cases 
by  row 

00 
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TABLE  XI 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  PRECIPITATION 
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TABLE  XII 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  CEILING 
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P  score  ean  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  seore  is  better. 

¥  scores  for  an  independent  data  sample  are  eomputed  using  probabilities  (relative  frequeney) 
obtained  from  eorresi>onding  developmental-sample  table. 


TABLE  XIV 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  VISIBIl  ITV 
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Q. 
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by  row 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better, 
scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 


TABLE  XV 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  VISIBILITY 


Total  eases 
by  row 
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.53 

.55 

o 

^.1 

1 

.32 

CO 

CO 

.45 

No.  of 

eases 

55 

588 

780 

1542 

901 

626 

134 

59 

4685 

VSBY  >  15 

13 

Relative 

frequency 

o 

00 

.59 

05 

lO 

CO 

.32 

.30 

o 

CO 

No.  of 

eases 

233 

1079 

971 

1032 

1  520 

374 

L  125 

24 

4358 

Y  range 
(lower  limit) 

o 

00 

Al 

o 

iq 

o 

CO 

- 1 

! 

o| 

.00  1 

-.20 

o 

in 

r 

i 

o 

in 

V 

Climatological 
totals  and 
frequencies 

ITS 


11 


03 

U 


I 

O 

u 

a 


C4 


a> 

u 

o 

I  a 


C/} 

> 


u  2 


0) 

K 


-H  I  a 


I 


58 


P  score  ean  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  seore  is  better. 

V  scores  for  an  independent  data  sample  are  eomputed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 


TABLE  XVII 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  TOTAL  CLD.  AMT. 


Total  cases 
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P  score  ean  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

V  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

“P  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 


TABLE  XIX 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  IOC 


Total  cases 
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69 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

V  scores  for  an  independent  data  sample  arc  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  PRECIPITATION 


Total  cases 
by  row 

147 

177 

553 

00 

00 

c- 

1770 

3159 

2962 

991 

10547 

OBSERVED  HSEOPITATICNCATEGORIES 

13 

Relative 

frequency 

No.  of 

cases 

^  1  13 

Relative 

frequency 

No.  of 

cases 

NO  PRECIP. 

13 

Relative 

frequency 

.53 

00 

00 

00 

c- 

o 

i  .  . 

.98 

.91 

No.  of 

cases 

09 

93 

397 

665 

1565 

2984 

o 

00 

00 

03 

971 

9615 

PRECIP. 

1—1  ^ 
a>  y 
>  ^ 

.C  a; 

1 1 

.59 

.47 

00 

CD 

Cvi 

.06 

^3 

OI 

o 

.09 

No.  of 

cases 

00 

84 

156 

123 

205 

ID 

t- 

00 

20 

932 

Y  range 
(lower  limit) 

in 

Al 

m 

CO 

CVJ 

00 

CD 

Oj 

o 

o 

o 

o 

V 

Climatological 
totals  and 
frequencies 

to 

It 

0) 


O 


c4 

O 

’S) 

O 


I 

o 

bc 

c 

o 

I  a 


64 


P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better, 
scores  for  an  independent  data  sample  arc  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSEHVED  VALUES  OF  CE'L'NG 
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TABLE  XXni 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  CEILING 
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TABLE  XXIV 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  VISIBILITY 
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© 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better 
P  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 


TABLE  XXV 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  VISIBILITY 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

V  scores  for  an  independent  data  sample  arc  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

V  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequeney) 
obtained  from  eorresponding  developmental-sample  table. 


TABLE  XXVII 

CONTINGENCY  TABLE  OF  INDEX  VALUES  vs.  OBSERVED  VALUES  OF  TOTAL 
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o 

to 

.20 
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.27 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

F  scores  for  an  independent  data  sample  arc  computed  using  probabilities  (relative  frequency) 
obtained  from  correspK)nding  developmental-sample  table. 


TABLE  XXVHl 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  IOC 
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For  definitions  of  IOC  Categories,  see  Table  XL. 

seore  (using  IOC  index  table  of  probabilities)  .56  p  score  (using  climatological  probabilities)  .60 


TABLE  XXIX 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  IOC 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

V  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  tabic. 


CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  CEILING 


Total  eases 
by  row 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

P  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental- sample  table. 
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TABLE  XXXIl 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  VISIBILITY 
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Total  cases 
by  row 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

V  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental -sample  table. 


TABLE  XXXIII 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  VISIBILITY 
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P  score  ean  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  seore  is  better. 

V  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 


TABLE  XXXIV 

CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  TOTAL  CLD.  AMT. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

V  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresp)onding  developmental -sample  table. 
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CONTINGENCY  TABLE  OF  INDEX  (Y)  VALUES  vs.  OBSERVED  VALUES  OF  TOTAL  CLD. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better, 
scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental- sample  table. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better, 
scores  for  an  independent  data  sample  arc  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental -sample  table. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

P  scores  for  an  independent  data  sample  arc  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  seorc  is  better, 
seorcs  for  an  independent  data  sample  are  eomputed  using  probabilities  (relative  frequeney) 
obtained  from  corresponding  developmental-sample  table. 
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P  score  can  range  from  a  minimum  of  0.0  to  a  maximum  of  2.0.  A  lower  score  is  better. 

P  scores  for  an  independent  data  sample  are  computed  using  probabilities  (relative  frequency) 
obtained  from  corresponding  developmental-sample  table. 
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